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Abstract — The  performance  of  HF  skywave  radar 
systems  is  customarily  referred  to  in  terms  of 
single-hop  propagation,  a  mechanism  which  pro¬ 
vides  illumination  of  the  earth’s  surface  out  to 
ranges  of  around  4000  kilometres.  In  practice, 
the  process  of  ionospheric  reflection  often  supports 
multiple  hops,  though  the  signals  are  inevitably 
subjected  to  much  greater  distortion  and  contam¬ 
ination.  In  this  paper  we  address  the  issue  of  ad¬ 
equacy  of  conventional  models  of  multi-hop  prop¬ 
agation.  We  formulate  a  detailed  model  which  ac¬ 
counts  for  intermediate  surface  scattering,  and  ob¬ 
tain  a  representation  in  terms  of  integrals  in  x  —  k 
space.  We  proceeed  to  evaluate  the  resulting  ex¬ 
pressions  for  several  cases  of  interest.  The  results 
demonstrate  that  it  is  vital  to  understand  the  com¬ 
plexities  of  multi-hop  propagation  if  this  method  of 
observation  is  to  be  exploited  for  remote  sensing  of 
the  ocean  at  extreme  ranges. 

Keywords — HF  radar,  remote  sensing,  ionosphere, 
radiowave  propagation 

L  Introduction 

HF  skywave  ’over-the-horizon’  radar  (OTHR)  ex¬ 
ploits  the  reflection  of  radiowaves  by  the  ionosphere 
to  illuminate  the  distant  earth’s  surface,  thereby  pro¬ 
viding  a  means  for  detecting  and  tracking  aircraft 
and  ships  at  ranges  up  to  several  thousand  kilome¬ 
tres.  In  addition,  information  about  the  sea  surface 
geometry  and  dynamics  is  embedded  in  the  radar 
echoes,  and  this  information  can  be  manipulated  to 
yield  estimates  of  a  variety  of  oceanographic  and  me¬ 
teorological  parameters  such  as  signiflcant  waveheight 
and  wind  velocity  [1].  Of  course,  the  vast  cover¬ 
age  achieved  by  OTHR  is  obtained  at  the  expense 
of  subjecting  the  signals  to  the  vagaries  of  the  iono¬ 
spheric  medium,  leading  to  distortion  and  contamina¬ 
tion  which  often  degrades  the  extracted  information 
[2]. 

The  zone  illuminated  by  radar  emissions  which  have 
undergone  a  single  reflection  from  the  ionosphere  ex¬ 
tends  typically  from  a  minimum  range  (’skip’  dis¬ 
tance)  of  about  1000  km  to  a  maximum  of  around 
4000  km,  and  for  most  purposes  this  constitutes  an 
upper  bound  to  the  range  of  observability  of  echoes  of 


Fig.  1.  Backscatter  ionogram  showing  one-,  two-  and 
three-hop  modes 

interest.  Nevertheless,  it  is  evident  that  signals  scat¬ 
tered  from  the  earth’s  surface  in  the  ’one-hop’  zone 
may  undergo  a  second  forward  ionospheric  reflection 
to  illuminate  more  distant  regions,  and  so  on. 

Figure  1  shows  a  backscatter  ionogram  in  which 
one-,  two-  and  three- hop  echoes  are  apparent.  Ray¬ 
tracing  codes  can  be  used  to  model  these  multi-hop 
paths,  as  illustrated  in  flgure  2,  which  was  generated 
with  a  two-dimensional  code  which  employs  the  modi- 
fled  Haselgrove  equations  to  compute  the  eikonal  rays 
through  a  model  ionosphere,  plotting  the  resulting 
trajectories  in  the  range- height  plane. 

Echoes  received  via  these  multi-hop  propagation 
paths  may,  in  principle,  be  analysed  to  provide 
oceanographic  information  out  to  ranges  in  excess  of 
6000  km,  noting  that  the  geometric  spreading  associ¬ 
ated  with  the  diverging  radar  beams  results  in  very 
coarse  spatial  resolution  at  long  ranges.  Of  course, 
the  increase  in  size  of  the  radar  resolution  cells  means 
that  the  clutter  power  increases  proportionately,  so, 
unlike  echoes  from  discrete  targets,  the  clutter  returns 
may  decay  relatively  slowly  with  range  and  hence  may 
maintain  the  clutter-to-noise  ratio  at  a  level  able  to 
support  analysis.  On  the  other  hand,  the  additional 
transits  through  the  ionosphere  cause  further  degra- 
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Fig.  2.  Ray  fan  computed  for  a  model  ionosphere,  showing 
multihop  modes 


Fig.  3.  Schematic  showing  the  progressive  diffusion  of  the 
transmitted  rays  on  successive  ground  reflections 


dation  of  the  signal,  so  the  availability  of  exploitable 
echoes  is  reduced. 

In  view  of  these  factors,  it  is  often  held  that,  un¬ 
der  clement  propagation  conditions,  the  exploitation 
of  multi-hop  echoes  can  proceed  in  a  manner  similar 
to  that  employed  for  one-hop  echoes,  using  essentially 
the  same  algorithms  for  parameter  extraction.  It  is 
the  purpose  of  this  paper  to  demonstrate  that  the  sit¬ 
uation  is  more  complex  than  appears  at  first  sight, 
and  that  analysis  of  multi-hop  echoes  demands  a  more 
sophisticated  model  of  the  radar  observation  process. 

II.  The  radar  process  model  for  multi-hop 

PROPAGATION 
A.  General  formulation 

The  process  model  equation  for  simple  one-hop 
propagation  can  be  written  [3] 


m  represents  internal  noise 

R  represents  the  receiving  complex,  including  anten¬ 
nas  and  receivers 
P  represents  the  signal  processing 
s  represents  the  signal  decomposition  after  processing 

The  general  form  of  the  extension  to  the  multi-hop 
propagation  case  can  be  written  as  a  sum  over  orders 
of  hop.  Immediately  we  observe  that  the  process  must 
account  for  scattering  from  areas  far  removed  from 
the  great  circle  path;  this  is  illustrated  schematically 
in  figure  3. 

It  will  be  convenient  in  what  follows  to  refer  on  oc¬ 
casion  to  the  number  of  ground  bounces,  rather  than 
the  number  of  ionospheric  hops,  particularly  to  facili¬ 
tate  description  of  non-reciprocal  paths.  The  process 
model  is  then  written  : 


s  —  pRMgSMj^Tw  +  pRM§n  -h  m  (1) 

where 

w  represents  the  selected  waveform 
T  represents  the  transmitting  complex,  including 
transmitters  and  antennas 

represents  propagation  from  transmitter  to  target 
scattering  region 

S  represents  all  scattering  processes  in  the  target  re¬ 
gion 

represents  propagation  from  target  scattering  re¬ 
gion  to  receiver 

n  represents  an  external  noise  source 

represents  propagation  from  noise  source  to  re¬ 
ceiver 


N  ns 

s  =  (2) 

nB=l  b=l 
N  rriB 

+  E  R[UMlitpSib)]M^^^^n  +  m 

nB='^  h=l 

Depending  on  the  spatial  filtering  properties  of  the 
operators  T,  M,  P  and  P,  some  of  these  terms  may  be 
suppressed.  Spatial  filtering  implemented  in  the  signal 
processing  operator  P  will  suppress  echoes  arriving 
from  directions  other  than  those  belonging  to  the  sub¬ 
space  defined  by  the  array  steering  vector.  Additional 
spatial  filtering  in  the  vertical  plane  arises  from  (i)  the 
elevation  radiation  patterns  of  the  antenna  elements, 
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embedded  in  the  transmit  and  receive  operators  T  and 
i?,  each  of  which  defines  an  intensity/sensitivity  dis¬ 
tribution  on  the  base  of  the  ionosphere,  and  (ii)  the 
propagation  operators  M,  which  incorporate  the  re¬ 
fractive  properties  of  the  ionosphere  and  hence  ac¬ 
count  for  focussing  and  defocussing.  Together,  these 
factors  usually  render  some  terms  in  the  expansion 
insignificant,  though  for  specific  applications  one  may 
need  to  verify  this. 

B.  Explicit  representation  of  the  operators 

In  order  to  proceed  with  calculations,  we  need  to 
choose  an  explicit  mathematical  representation  for  the 
operators.  The  choice  depends  on  the  physical  domain 
in  which  the  radar  process  is  to  be  studied.  For  exam¬ 
ple,  if  the  transformation  of  the  polarisation  state  of 
the  signal  were  the  only  feature  of  interest,  the  anal¬ 
ysis  could  be  carried  out  using  the  Jones  calculus. 

It  is  natural  and  convenient  to  adopt  a  representa¬ 
tion  based  on  integrals  over  the  space  and  wavevector 
variables;  for  the  present  we  shall  suppress  the  (evolv¬ 
ing)  polarisation  state  of  the  field  and  the  temporal 
dependence  of  the  operators. 

We  can  then  write,  to  second  order  : 


f  =  —PTGT{k){.) 

(3) 

II 

^1 

yr; 

3 

(4) 

S  =  J  dk5{k' -k'.h-k  +  k.h)5{r.h  +  k.h){.){5) 

+  J  dKlF^^\ki,ks,Kl)S{Kl){.) 

+  J  J  dKldK2F^^\ki,ks,Kl,K2)S{Kl)S{K2){-) 

R^f^GR{k){.)  (6) 

where  G{k' ,r^ ,k,r)  is  the  Green’s  function  for  sky- 
wave  propagation  and  S{n)  denotes  the  directional 
wave  spectrum  of  the  sea  surface.  F^{ki^ks^  ni)  and 
F‘^{ki^ks^  1^2)  denote  the  kernels  for  the  first  and 
second  order  scattering  respectively. 

C.  Construction  of  the  model  for  two-hop  propagation 

The  generality  of  the  formal  model  enables  it  to 
incorporate  diverse  effects  including  dispersion,  gen¬ 
eralised  Faraday  rotation,  tilts  and  ionospheric  distur¬ 
bances  and  so  on.  Our  focus  here  is  on  the  relative  im¬ 
portance  of  various  contributions  to  the  received  field. 
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Fig.  4.  Theoretical  Doppler  spectra  of  sea  clutter  for  var¬ 
ious  categories  of  scattering  geometry 


so  we  shall  simplify  the  physics  as  much  as  possible 
without  compromising  the  essential  aspects.  This  has 
the  advantage  of  speeding  up  the  calculations  -  on  a 
2GHz  Pentium  PC,  the  model  can  track  400,000  rays 
in  about  10  minutes. 

For  the  computations  reported  here,  we  adopt  a 
spherical  earth,  a  concentric  spherical  mirror  as  our 
ionosphere,  consider  only  the  zeroth  and  first  order 
terms  in  the  expansion  for  the  sea  scatter,  assume 
a  land-free  ocean,  and  choose  a  transmit  antenna 
pattern  with  uniform  sector  gain  in  both  the  verti¬ 
cal  and  horizontal  planes.  Receive  beams  are  formed 
over  360°  assuming  element  patterns  isotropic  in  az¬ 
imuth  with  uniform  sector  gain  in  elevation.  Array 
beamwidth  is,  naturally,  a  function  of  azimuth.  The 
sea  surface  is  assumed  to  be  fully-developed  over  the 
frequency  band  corresponding  to  the  Bragg  resonant 
scattering  mechanism,  and  isotropic  in  azimuth. 

Ultimately  the  distribution  of  energy  in  the  Doppler 
domain  is  of  central  interest,  so  it  is  important  to 
model  the  scattering  from  the  sea  surface  with  ade¬ 
quate  fidelity.  For  this  purpose  we  employ  the  JIN- 
SCAT  code  developed  in  DSTO  for  remote  sensing 
and  ship  detection  applications.  Representative  spec¬ 
tra  computed  with  JINS  CAT  for  several  broad  classes 
of  scattering  geometry  are  shown  in  figure  4. 

While  the  full  Doppler  spectrum  can  be  computed 
in  an  accurate  and  natural  way  with  our  model,  the 
present  paper  restricts  attention  to  the  total  power, 
without  a  decomposition  in  Doppler.  As  shown  below, 
it  turns  out  that  computing  just  the  power  suffices  to 
answer  the  key  questions  posed  here;  results  generated 
for  the  higher-dimensional  analysis  will  be  published 
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elsewhere.  A  further  restriction  imposed  here  is  to  as¬ 
sume  collocated  transmit  and  receive  antennas;  this  is 
done  simply  to  facilitate  comparison  with  common  ex¬ 
perience  -  the  model  is  equally  applicable  to  arbitrary 
bistatic  geometries  [4]. 

We  shall  consider  the  case  of  a  radar  employing  two 
uniform  linear  arrays,  a  transmitting  array  10  A  in 
length  and  a  receiving  array  lOOA  in  length.  For  a 
specified  ionospheric  height  h,  and  a  take-off  elevation 
angle  9  from  point  r,  the  ground  and  slant  ranges  to 
the  ground  bounce  point  r'  are  given  by  : 


a  — 

Re  *  5m(|  -h  9) 

(7) 

Re  T  h 

p  = 

—  —  9  —  arcsin((T) 

(8) 

Rgr  — 

2  *  i?e  *  /? 

(9) 

Rsi  = 

_  sind 

2  *  Re  * 

a 

(10) 

where  Re  is  the  radius  of  the  earth  .  It  is  convenient 
to  employ  a  mixed  coordinate  system,  wherein  k  is 
defined  relative  to  the  local  normal  h  at  r,  while  k' 
is  defined  relative  to  the  local  normal  ft'  at  This 
facilitates  the  calculation  when  the  surface  scattering 
is  the  primary  object  of  study. 

The  Green’s  function  then  takes  the  form 

G{k\r  ^k^f)  =  S(k' —k)6{k'  .n  -\-k.h)  (11) 

Rsi 

We  note  that  most  of  the  process  model  operators 
depend  on  elevation  angle,  including  not  only  the  an¬ 
tenna  patterns  but  also  the  scattering  from  the  ter¬ 
restrial  surface.  As  shown  in  Figure  5,  the  elevation 
angles  of  one-and  two-hop  modes  can  be  very  differ¬ 
ent  at  a  given  range,  which  can  lead  to  substantial 
changes  in  the  scattering  coefficient  when  the  Bragg 
frequency  is  close  to  the  wave  spectrum  peak. 

III.  Coordinate  registration  for  multi-hop 

PROPAGATION 

The  visual  appearance  of  the  backscatter  ionogram 
presented  in  figure  6  suggests  that  the  two-hop  mode 
extends  the  range  coverage  from  the  one-hop  limit  at 
3500  km  to  a  new  maximum  of  about  5500  km.  This 
is  not  strictly  correct  -  the  Y-axis  in  the  figures  1 
and  6  is  the  ’group’  or  ’slant’  range,  and  the  conver¬ 
sion  to  ground  range  must  take  account  of  the  actual 
propagation  path.  In  the  case  of  two-hop  modes,  the 
additional  distance  travelled  in  the  course  of  under¬ 
going  two  ionospheric  reflections  must  be  taken  into 


Elevation  angles  of  arrival  for  1-hop  and  2-hop  modes 


Fig.  5.  Variation  of  elevation  angle  with  ground  range  for 
one-  and  two- hop  modes 


Fig.  6.  Backscatter  ionogram  showing  the  transition  from 
the  one-hop  regime  to  the  two-hop  regime 


account.  The  relevant  relationships  and  corrections 
are  presented  in  figures  7-9.  Figure  7  plots  the 
slant  (group)  range  as  a  function  of  ground  range, 
parametrised  by  ionospheric  height,  for  both  one-  and 
two-hop  propagation.  This  information  can  be  used 
to  calculate  the  two-hop  range  correction  factor,  as 
presented  in  Figure  8,  which  shows  the  range  differ¬ 
ence  of  two-hop  as  opposed  to  one-hop  propagation  to 
a  given  range.  Here  and  in  all  figures,  the  horizon  for 
a  given  ray  is  indicated  by  termination  of  the  corre¬ 
sponding  curve.  Finally,  it  is  useful  to  plot  one  mode 
against  the  other,  as  is  done  in  Figure  9,  to  empha¬ 
sise  the  relative  (as  opposed  to  the  absolute)  range 
discrepancy. 
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Group  range  to  a  common  ground  point  for  one-hop  and  two-hop  modes 


Fig.  7.  Comparison  of  group  ranges  for  one- hop  and  two- 
hop  propagation  to  a  given  ground  range 


Difference  between  1-hop  and  2-hop  group  ranges  to  a  common  ground  range 


Fig.  8.  The  range  difference  introduced  by  two-hop  prop¬ 
agation  as  opposed  to  one-hop  propagation 


IV.  Climatology  of  multi-hop  propagation 

The  Jindalee  radar  frequency  management  system 
records  backscatter  ionograms  routinely,  twenty  four 
hours  a  day,  and  has  done  so  since  1982  with  some 
short  interruptions.  Until  1999  the  data  was  acquired 
one  beam  at  a  time  (there  are  eight  beams  spanning 
the  90  deg  arc  to  the  north-west  of  Australia)  but  since 
July  1999,  when  an  eight-channel  receiver  was  fitted, 
all  beams  have  been  collected  simultaneously.  Typi¬ 
cally  the  full  set  of  ionograms  is  updated  every  4  min¬ 
utes,  covering  ranges  out  to  12,000  kilometres. 

This  extensive  database  is  ideally  suited  to  the  sta¬ 
tistical  analysis  of  the  directional,  diurnal,  seasonal 


2-hop  versus  1-hop  group  range  to  a  common  ground  range 


Fig.  9.  Group  ranges  compared  :  two-hop  versus  one-hop 


and  sunspot  number  dependence  of  multihop  propa¬ 
gation.  Such  an  analysis  lies  beyond  the  scope  of  the 
present  paper,  but  we  note  in  passing  that  in  prepar¬ 
ing  this  paper  we  examined  a  single  month’s  data  in 
some  detail  and  identified  distinct  patterns  in  multi¬ 
hop  coverage. 

V.  Results 

The  first  property  of  interest  is  the  extent  of  the 
spread  in  azimuth  of  the  radar  echoes  after  diffuse 
scattering  from  the  first  bounce  footprint.  Figure  10 
shows  a  histogram  of  azimuthal  angle-of-arrival  of  10^ 
rays  originally  transmitted  in  a  broad  beam  of  full 
width  120°  centred  on  270°.  It  is  apparent  that  the 
tendency  is  for  considerable  broadening  but  with  the 
distribution  centred  around  the  axis  of  the  transmit¬ 
ted  beam.  This  is  a  consequence  of  the  fact  that  the 
diffusely-scattered  rays  are  generated  with  a  uniform 
distribution  in  azimuth  from  the  first  bounce  point 
onwards.  Note  that  this  figure  reflects  only  the  angle- 
of-arrival  of  the  rays,  not  their  amplitudes. 

The  key  issue  is  the  magnitude  of  the  contribution 
from  the  diffusely-scattered  field  relative  to  that  of 
the  specularly  scattered  field.  It  is  appropriate  to  use 
a  common  scale  when  presenting  range-Doppler  maps 
of  the  synthesised  receiver  output,  and  this  was  first 
implemented  by  using  the  global  maximum  over  the 
families  of  specularly-reflected  and  diffusely-scattered 
triple-bounce  modes.  Later  an  alternative  approach 
was  adopted  -  we  included  in  the  output  the  returns 
arriving  via  singly-scattered  one-hop  propagation,  and 
used  the  global  maximum  over  all  three  families  as 
the  scaling  constant.  This  approach  facilitates  a  com- 


501 


7000 


Fig.  10.  A  histogram  of  azimuthal  angle-of-arrival  at  the 
receiving  array  for  triple-bounce  diffusely- scattered 
multihop  rays 


parison  between  one-hop  and  two-hop  specularly  re¬ 
flected  echoes.  An  example  of  a  ’triptych’  display 
of  the  range-Doppler  maps  for  the  three  families  of 
rays  is  shown  in  Figure  11.  For  ease  of  interpretation, 
only  10  transmitted  rays  are  shown,  with  each  yield¬ 
ing  a  single  received  ray  for  the  non-diffusely  scattered 
modes.  In  contrast,  the  diffusely-scattered  rays  spread 
more  or  less  uniformly  over  a  substantial  range  extent. 

The  relative  magnitude  of  the  diffusely-scattered  re¬ 
ceived  field  relative  to  the  specularly  reflected  two-hop 
return  is  seen  to  be  around  -30  dB  for  the  assumed 
parameter  values.  This  contribution  would  fall  be¬ 
low  the  subdominant  Bragg  line  on  almost  all  real¬ 
isations  of  the  sea  clutter  Doppler  spectrum,  but  is 
certainly  comparable  with  much  of  the  second-order 
clutter  spectrum.  It  is  the  latter  which  is  central  to 
the  extraction  of  most  of  the  detailed  sea  state  infor¬ 
mation  sought  in  HF  radar  remote  sensing. 

VI.  Conclusion 

We  have  developed  a  computational  model  able  to 
compute  the  contributions  of  diffusely-scattered  fields 
arising  from  non-specular  reflection  (bistatic  scatter¬ 
ing)  at  the  ground  reflection  points  for  multiple-hop 
propagation  paths.  The  significance  of  these  echoes  is 
determined  by  their  magnitude  relative  to  the  echoes 
received  via  specular  reflection  at  the  bounce  points, 
and  by  their  spectral  content. 

Early  calculations  reveal  that  the  diffusely- 
scattered  returns  lie  SOdS  below  the  specularly- 
reflected  returns  for  the  two-hop  propagation  case. 


This  results  has  important  implications  for  remote 
sensing  : 

(a)  under  conditions  corresponding  to  the  assumed 
parameter  values,  the  diffusely-scattered  returns  are 
not  so  strong  as  to  bias  the  Bragg  line  ratio  apprecia¬ 
bly,  so  information  such  as  wind  direction  estimated 
from  the  Bragg  line  ratio  should  not  be  affected 

(b)  procedures  for  estimating  sea  state  and  other 
ocean  surface  characteristics,  including  inversion  to 
yield  the  full  directional  wave  spectrum,  may  be  seri¬ 
ously  corrupted  by  the  presence  of  diffuse  multi-hop 
contributions 

(c)  the  relative  magnitude  of  the  diffusely-scattered 
contributions  is  strongly  influenced  by  the  antenna 
patterns,  both  transmit  and  receive,  with  obvious  im¬ 
plications  for  radar  design 

(d)  the  contributions  from  different  areas  will,  in 
general,  have  experienced  different  phase  path  mod¬ 
ulations  as  a  consequence  of  the  spatial  and  tempo¬ 
ral  variability  of  the  ionosphere,  leading  to  additional 
contamination  of  the  desired  Doppler  spectra 

(e)  the  prospect  of  reliable  remote  sensing  at  ex¬ 
treme  ranges  via  multi-hop  propagation  will  depend 
on  the  development  of  new  techniques  able  to  de¬ 
tect,  estimate  and  perhaps  even  compensate  for  dif¬ 
fuse  scattering  contamination. 
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Range-Doppler  map  for  (a)diffuse,  (b)one-hop  and  (c)two-hop  propagation 
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Fig.  11.  Range-Doppler 
(backscatter),  and  (c) 


maps  computed  for  :  (a)  diffusely-scattered  triple-bounce  rays,  (b)  singly-scattered  rays 
specularly-scattered  triple-bounce  (two-hop)  rays 
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